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GRAPHICAL SOLUTION OF A DIFFERENTIAL EQUATION. 279 

If we have (D — a) 2 y = 0, we may therefore get (D — a)y = Cie ax . On 
repeating this operation, 

y = e ax Je" w cie° x dx = c x xe ax + c e ax . 

If we have (D — a) s y = 0, we get in precisely the same way 

y = C2X 2 e ax + Cixe ax + c e a *. 

The complementary function would then be completely known on proving the 
usual theorem that if v u v 2 , ■• • v n are solutions of f(D)y = then ci»i + czv 2 
+ • • • c„v„ is also a solution. 

One need hardly draw attention to the fact that the complete value of y will 
be obtained if the factors of 

(D - ai)(D -a 2 ) ■■• (D- a n )y = X 

are removed in turn by the above rule, the constant of integration being intro- 
duced each time. 



A GRAPHICAL SOLUTION OF THE DIFFERENTIAL EQUATION 

OF THE FIRST ORDER.* 

By T. R. RUNNING, University of Michigan. 

The following graphical solution of the differential equation of the first order 
was derived while trying to solve a problem involving two such equations. As 
an illustration, these equations will be solved after a discussion of the method of 
solution has been given. 

If the primitive of a differential equation which expresses the condition of 
a problem can not be found, or can be found with difficulty, it may i e desirable 
to have a method of approximate solution which does not involve a great deal 
of labor. It is thought that the method given is more easily carried out than 
any other which has come to the writer's attention. 1 

The solution depends upon the simple fact that the area under the derived 
curve is equal to the corresponding ordinate of the integral curve. In Fig. 1 

Area OPQ = RQ. 
If in the equation 

we assign particular values to y and plot the resulting equations with dy/dx as 

* First part of a paper presented to the American Mathematical Society at Madison, Sep- 
tember 9, 1913. 

1 For a comparison of approximate methods see " Beitrag zur naherungsweisen Integration 
totaler Differentialgleichungen," by W. Kutta, Zeitschrift ftir Mathematik und Physik, Vol. 46, 
pp. 435-453. 
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ordinate and x as abscissa, we shall have, designating each curve by the corre- 
sponding value of y, a family of curves as shown in Fig. 2. 
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Fig. 1. 
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Fig. 2. 



Suppose the initial condition to be y = y 2 when x = 0. Draw a curve, 
which for an approximation may be a straight line, from P , where the curve y 2 
crosses the axis of dy/dx, to the curve y s , such that the area under it shall be 
equal to y 3 — y%. From this point, Pi, draw a line to a point P 2 on the curve 
2/ 4 such that the area under the line PiP 2 shall be equal to t/ 4 — yz- This process 
continued will give any desired number of points on the integral curve. The 
curve connecting the points P Pi • • • P» is the derived curve. 

It is to be noted that the integral curve is not obtained by graphically inte- 
grating the derived curve, but the coordinates of its points are read from the 
figure. Thus (x , 2/2), («i, ys), (%2, yt)> ( x s, y&), etc., are points which determine 
the integral curve. If a close approximation is required the values of y must 
be so chosen that the corresponding curves will lie near together. 

The problem mentioned was one which gave rise to the two differential 
equations 



(1) 



(2) 



dy 
dx 



dy _ 1800a; - 132y* 
dx ~ lit/* 4- 1000 ' 

1740 - 240a- - 132y» 



\ly* + 1000 



It was required to find the maximum value of y and the corresponding value 

of x, knowing that when x = 0, y = 0, equation (1) was to be satisfied between 

29 29 

x = and x = ^r, equation (2) between x = ^7 and x = 7j. There was another 

equation to be satisfied beyond x = 7J but, as the conditions clearly showed, 
the maximum value of y occurs before x reaches the value 7J, it need not be 
considered. 

The curves which correspond to different values of y are straight lines and 
shown in Fig. 3. 
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Fig. 3. 

y will be a maximum when dyldx = 0. This value is 2.926 and the corre- 
sponding value of x is 4.492. The values chosen for y are indicated on the 

29 
figure. It is seen that when x = 57, y = .623. This value must, therefore, be 

used in both differential equations. 

The method was applied to the problem given in Kutta's article, above 
referred to, for comparing the five different approximations. With squared 
paper ruled to twentieths of an inch and using two inches as the unit, the approxi- 
mation was closer than that of any one of the first four but not quite so close as 
Kutta's, while the work necessary to obtain the result was much less. 

In a later paper it is proposed to extend the method to the solution of the 
differential equation of the second order. 
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W. H. Bussey, Chairman of the Committee. 

Practical Mathematics. By Claude Irwin Palmer, Assistant Professor of 

Mathematics in the Armour Institute of Technology. The McGraw-Hill 

Book Company, New York, 1913. Four volumes. $0.75 each. 

We have here in four small, neatly bound and clean cut volumes a very 
interesting presentation of elementary mathematical subjects from Arithmetic 
through Trigonometry. The work has been prepared for persons of maturity, 
being, as stated by the author, the outgrowth of a course in practical mathematics 
given in the evening classes of the Armour Institute of Technology. 

The titles of the four little volumes are: 

Part I. Arithmetic with Applications (150 pages). 

Part II. Geometry with Applications (166 pages). 

Part III. Algebra with Applications (189 pages). 



